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Abstract. We present &BVI photometric and astrometric catalogue of the open cluster NGC 3960, down to limiting magnitude 
V ~ 22, obtained from observations taken with the Wide Field Imager camera at the MPG/ESO 2.2 m Telescope at La Silla. 
The photometry of all the stars detected in our field of view has been used to estimate a map of the strong differential reddening 
affecting this area. Our results indicate that, within the region where the cluster dominates, the E(V — I) values range from 
0.21 up to 0.78, with E(V - 1) = 0.36 (E(B - V) = 0.29) at the nominal cluster centroid position; color excesses E(V - I) 
up to 1 mag have been measured in the external regions of the field of view where field stars dominate. The reddening corrected 
color-magnitude diagram (CMD) allows us to conclude that the cluster has an age between 0.9 and 1.4 Gyr and a distance 
modulus of {V - M v )o = 11-35. 

In order to minimize field star contamination, their number has been statistically subtracted based on the surface density map. 
The empirical cluster main sequence has been recovered in the V vs. V — I and in the J vs. J — Ks planes, using optical 
and infrared data, respectively. From these empirical cluster main sequences, two samples of candidate cluster members were 
derived in order to obtain the luminosity distributions as a function of the V and J magnitudes. The Luminosity Functions have 
been transformed into the corresponding Mass Functions; for M > 1 M©, the two distributions have been fitted with a power 
law of index av = 2.95 ± 0.53 and aj — 2.81 ± 0.84 in V and in J, respectively, while the Salpeter Mass Function in this 
notation has index a — 2.35. 

Key words. Open clusters - individual: NGC 3960 - photometry - astrometry - differential reddening - Luminosity and Mass 
Function 



. 1. Introduction 

Old Galactic open clusters are stellar systems which can 
survive only if they have a large initial populat ion and lie 
far from dense molecular clouds Jvan den Bergh & McCluret 
1980: lde La Fuente Marcosl! 19971 For this reason they are less 
numerous than the young open clusters and are mainly found in 
the external region of the Galaxy, where there is a lower likeli- 
hood of catastrophic encounters. 

Cluster stars have the same distance and chemical compo- 
sition; for this reason stellar clusters are among the most ap- 
propriate objects to study the Initial Mass Function (IMF), that 
is one of the most crucial ingredients for dynamical evolution 
models and for galaxy formation and stellar evolution models. 
Comparison of the IMF of old open clusters with the IMF of 
nearby and young open clusters is fundamental to look for sim- 
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ilarities or discrepancies of the star formation processes in the 
Galaxy disk, in different epochs and in different environments. 

Accurate determination of the open cluster IMF is, how- 
ever, a difficult task, because these objects are in general poorly 
populated and, due to their location in the thin disk, are highly 
contaminated by field stars seen along the same line of sight. 
In the case of old open clusters, these difficulties are increased, 
as they are very distant from the Sun and thus often strongly 
obscured due to interstellar absorption. High-quality photome- 
try and precise extinction determinations are crucial to obtain a 
correct estimate of the IMF. 

With an age of about 1.2 Gyr and a small spatial concentra- 
tion, NGC 3960 is an interesting example of an old open clus- 
ter suitable for IMF studies. It is a low-latitude cluster, located 
at the celestial and Galactic coordinates RA=1 l' l 50 m 54 ;5 and 
Dec=-55°42'(J2000) and I =294°42 and b =+6°. 17, respec- 
tively, about 1850 pc from the Sun. 

Ph ysical parameters of NGC 3960 were estimated bv ljanesl 
dl98ll) . through photoelectric BV and intermediate-band David 
Dunlop Observatory (DDO) photometry. One blue and one vi- 
sual photographs were obtained in order to extend the mag- 
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nitude sequences to fainter limits. From the BV photometry, 
the author concludes that the cluster is (0.5 — 1.0) x 10 9 yr 
old, while from the DDO photometry for several cluster giants 
and using the Janes] i 19771) method, he derives the reddening 
E(B -V) = 0.29 ± 0.02. A distance modulus {V - M v ) = 
11.1 ± 0.2 and metallicity [Fe/H]=— 0.30 have also been esti- 
mated by the same author. [Fe/H l=— 0. 34 ± 0.08 was subse- 
quently derived by iFriel & Janesl Jl993l) from a spectroscopic 
study of 7 giant stars of NGC 3960. 

Because of its low-latit ude position , in the direction of 
the Carina spiral feature (Miller, 119721) . effects of differen- 
tial reddening are expected across the region of t he sky in- 
cluding the c l uster members. Feitzin ger & Stuewel if 1984) and 
IHartlev et alJ Jl986l) indicate the presence in their catalogues 
of two dark nebulae f alling i n the field of NGC 3960. More re- 
cently, iDufrT&Bici (|2002) gave values E(B - V) = 0.51 
and E(B — V) = 0.59 for these two overlapping dark nebu- 
lae, respectively, derived from a 100 ^m dust thermal emission 
map. These regions are located about 13 arcmin from th e clus- 
ter c enter, where the E(B — V) value is about 0.29 ( I Janesl 
Il98lh . This non negligible difference in the E(B — V) value, 
therefore, suggests the existence of variable extinction across 
the field of NGC 3960. 

In this work, usin g a procedure similar to t hat used by 
IPiotto et all i 19991) and lvon Braun & Mateol l l200ll) . we attempt 
to construct an extinction map of our field of view, in order to 
correct the photometric catalogue of candidate cluster members 
for differential reddening effects. 

In Section|2] we describe the observations and the data re- 
duction procedure used to obtain the photometric/astrometric 
catalogue; in Section [5] we present the CMD and the method 
used to obtain the extinction map and the reddening corrected 
catalogue. In Section|4]we describe the method adopted to de- 
rive the cluster parameters, while in Section[5]we describe the 
statistical procedure used to recover the empirical cluster main 
sequence from optical and infrared data and the photometric se- 
lection of candidate cluster members. In Section[6] we present 
the Luminosity Functions obtained from optical and infrared 
data and the corresponding Mass Functions. Finally, we sum- 
marize and discuss our results in Sectional 

2. Cluster BVI Photometry and Astrometry 

2. 1. Observations and Data Reduction 

NGC 3960 is an open cluster observed in the Pre-FLAMES 
Survey (PFS), carried out as part of the ESO Imaging Survey 
(EIS), aimed to provide BVI imaging data for use in connection 
with FLAMES, a multi-fiber spectrograph instrument on the 
ESO VLT UT2 (KUEYEN) Telescope. The observations for 
this survey were collected using the Wide Field Imager (WFI) 
camera mounted at the Cassegrain focus of the MPG/ESO 
2.2 m Telescope at La Silla (Chile). This instrument consists of 
a 4 x 2 mosaic of 2k x 4k CCD detectors with narrow inter-chip 
gaps of width 23"8 and 14"3 along right ascension and decli- 
nation, respectively, yielding a filling factor of 95.9%. With a 
pixel size of 0.238", each chip covers a field of view of 8' 12 x 
16(25, while the full field of view is 34' x 33'. 



According to the PFS observing strategy, the images were 
obtained in the BVI pass-bands; for each filter, one short ex- 
posure of 30 seconds, useful to avoid saturating bright objects, 
and two deep exposures, of 4 minutes each, have been taken. 
In order to cover the inter-chip gaps, these observations were 
dithered by 30" both in right ascension and declination. Table 
Ogives the log-book of the observations and Figure^shows a 
deep image of our field of view in the V band. The cluster is 
approximately located in the central part of the field of view 
where it is barely visible within a radius of about 7 arcmin. 
Strong effects of differential reddening cause a non-uniform 
spatial distribution of the field stars in the external regions of 
the field of view. In addition to the science target exposures, 
a set of technical frames for the instrumental calibration were 
obtained during the same night. 

To complement the optical data, we have utilized JHKs 
infrared (IR) photometry from the All-Sky Point Source 
Catalogue of the Two Micron All Sky Survey (2MASS) 1 
JCarpenteilEoOll) available on the WEB . 

The instrumental calibration of the WFI CCD images has 
been performed using the mscred package, a mosaic specific 
task implemented as an IRAF package for the NOAO Mosaic 
Data Handling System (MDHS). 

The first stage of the data reduction process has been the 
removal of the instrumental signatures from each raw CCD im- 
age. First of all, we have subtracted the electronic bias using the 
overscan region that has been subsequently removed. Exposure 
bias patterns have been subtracted using a zero calibration mo- 
saic obtained as the median of two bias frames. Finally, fiat 
fielding for each filter has been performed using a set of dome 
flat fields combined into a median master flat. Due to some sat- 
urated points, we could not use the sky flat fields. 

Images in the / band have required a special treatment be- 
cause of the strong effects of fringing. In order to remove this 
instrumental artifact we have subtracted the fringing pattern 
provided by the MPG/ESO 2.2 m Telescope team 3 , scaling it 
to each exposures. 

Instrumental magnitudes were obtained for each chip of 
each mosaic im age by using the DAOPHOT II /ALLSTAR 
JStetsoni Il987h and ALLFRAME JStetsorl 1 1994 photomet- 
ric routines. The point-spread functions (PSFs), determined for 
each chip and filter, take into account the variation in the PSF 
across the field of view of each WFI chip. 

In order to obtain the aperture correction to the profile- 
fitting photometry, growth curves were determined. First, 
we have selected the stars used to define the PSF model; 
next, all other objects were removed from the frames and 
aperture phot ometry was c arried out at a variety of radii. 
DAOGROW (Stetson, 1990) was used to derive growth curves 
and COLLECT to calculate the "aperture correction" coeffi- 
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Table 1. Log-book of the observations. 



Target 


RA (J2000) 


Dec (J2000) 


Night Filter 


Exp. Time 


seeing 


EIS name 


(h m s) 


(d m s) 




[s] 


FWHM["] 


NGC 3960 


11 50 55.0 


-55 41 36 


24-25 Feb 2000 B 


1 x 30 + 2 x 240 


0.99-1.21 


(OC21) 






V 


1 x 30 + 2 x 240 


0.86-1.38 








I 


1 x 30 + 2 x 240 


0.81-1.19 




N ■ 



Fig. 1. A deep WFI image of the region around NGC 3960 taken in the B filter. The field center is RA=ll ,l 50 m 51 f :6 and 
Dec=-55°42'07"(J2000). Each chip is 8(12 x 16(25, while the full field of view is 34' x 33'. The cluster is approximately 
located in the central part of the field of view where it is barely visible within a radius of about 7 arcmin. Strong effects of 
differential reddening cause a non-uniform spatial distribution of the field stars in the external regions of the field of view. 



cient for each chip from the difference between PSF and aper- 
ture magnitudes of the selected stars. 



Calibration to the Johnson-Cousin photometric system was 
performed by means of a set of lLandoltllll992l) standard fields 
(SA98, SA104, PG0918 and Ru 152), observed during the 
same night through the 8 chips. 



In order to derive the transformation coefficients to the stan- 
dard system we have used equations of the form: 

v = V + A Q + A 1 x X + A 2 x (B-V), 

b =B + B + B 1 xX + B 2 x {B-V), (1) 

i =I + C + C\xX + C 2 x(V-I). 

where v, b and i are the standard star instrumental magnitudes, 
corrected to the exposure time of 1 s; X is the airmass and V, 



4 



Prisinzano et al.: NGC 3960 



Table 2. Coefficients of the transformation to the standard system for each filter and for each chip. Col. 3 gives the number of 
standard stars used to derive the zero points and the color terms (Col. 4 and 5), imposing the mean extinction coefficients of La 
Silla (Col. 6); finally, Col. 7 gives the corresponding average standard errors. 



Filter 


Chip# 


N. stars. 


Zero Point 


Color Term 


Extinction 


Av. St. Error 


V 


1 


24 


0.987 ± 0.007 


-0.097 


0.14 


0.034 




2 


306 


1.022 ±0.002 


-0.097 


0.14 


0.031 




3 


136 


1.062 ± 0.003 


-0.097 


0.14 


0.030 




4 


4 




-0.097 


0.14 






5 







-0.097 


0.14 






6 


51 


1.095 ± 0.005 


-0.097 


0.14 


0.033 




7 


67 


1.048 ± 0.003 


-0.097 


0.14 


0.022 




8 


11 


0.982 ± 0.007 


-0.097 


0.14 


0.024 


B 


1 


25 


0.668 ± 0.012 


0.249 


0.25 


0.053 




2 


328 


0.723 ± 0.002 


0.249 


0.25 


0.038 




3 


148 


0.772 ± 0.002 


0.249 


0.25 


0.028 




4 


4 




0.249 


0.25 






5 







0.249 


0.25 






6 


49 


0.810 ±0.008 


0.249 


0.25 


0.054 




7 


56 


0.754 ± 0.004 


0.249 


0.25 


0.030 




8 


11 


0.697 ±0.019 


0.249 


0.25 


0.062 


I 


1 


22 


1.989 ±0.012 


0.130 


0.09 


0.056 




2 


248 


2.042 ± 0.002 


0.130 


0.09 


0.035 




3 


123 


2.063 ± 0.003 


0.130 


0.09 


0.038 




4 


6 




0.130 


0.09 






5 







0.130 


0.09 






6 


42 


2.072 ± 0.007 


0.130 


0.09 


0.047 




7 


47 


2.066 ± 0.005 


0.130 


0.09 


0.033 




8 


13 


2.009 ± 0.007 


0.130 


0.09 


0.026 



B and / are the magnitudes in the standard system. We first 
attempted to derive all the coefficients (the zero points, the ex- 
tinction and the color terms) for each chip, but by analyzing 
the photometric residuals as a function of the time through the 
night, we noticed that when the PG0918 standard field was ob- 
served, the night was not photometric, as the mean residuals 
systematically change with time. Since the problem does not 
affect the other science observations, we have discarded all the 
observations of the PG0918 field and we have only used the 
Landolt fields SA98, SA104 and Rubin 152. In order to im- 
prove the photometric calibration, the number of the standard 
stars in each of thes e fields has bee n increased using secondary 
standards defined in Stetson (2000). Unfortunately, neither the 
Landolt's nor the Stetson's catalogue cover the chips 4 and 5 of 
our field with a sufficient number of standards to determine in- 
dependent coefficients for these chips. Neighboring chips have 
been used to calibrate science observations in these chips, as 
described at the end of this section. Imposing mean extinction 
coefficients (Stetson, private communication), a zero point and 
a color term were determined separately for the chips where a 
sufficient number of standard observations were available. We 
have found that the color terms computed for the chips 2, 3, 6 
and 7 are consistent with each other, with variations < 0.01 in 
B, < 0.02 in V and < 0.05 in I; for the chips 1 and 8 we have 
found different values because they have been computed on a 
smaller number of standard stars (see Table|2j. Considering the 
homogeneous characteristics of the instrument, we have chosen 



to compute new zero points fixing the color term to the median 
value of the color terms of the chips for which the number of 
measured standard stars was sufficiently large. The resulting 
coefficients are given in Table |2] together with the correspond- 
ing average standard errors. Points above the 3 a level were 
not considered for the calibration. We note that the variations 
among the photometric zero points of different chips (Table|2j, 
are < 0.11 in V, < 0.1 4 in B and < 0.08 i n I, in agreement 
with the results found in IZoccali etalJj2003h . 

As already mentioned, the NGC 3960 field observations 
were dithered both in right ascension and declination to cover 
the inter-chip gaps. This means that contiguous chips have sev- 
eral stars in common. We have compared the magnitudes of 
these stars to verify consistency of our photometry among dif- 
ferent chips. The mean value and the root mean square of the 
residuals for the stars with V < 20 are given in Table. [3] 
Finally, the photometric zero points of the chips 3 and 6 have 
been used to calibrate the stars in the chips 4 and 5, respec- 
tively. Corrections < 0.02, < 0.01 and < 0.02 mag to these 
photometric zero points have been found, for the magnitudes 
V, B and /, respectively, using the photometric comparison of 
the common stars in the overlapping regions. 

Completeness of our star list and accuracy of our photom- 
etry for each chip have been determined by tests with artificial 
st ars; a total of about 13 50 stars have been added as described 
in Prisinz aiio et alJ ll2003l) . Our data are 100% complete above 
\ ' = 20 and more than 50% complete above V ~ 21, 1 ~ 19.5 
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Table 3. The mean value and the root mean square of the V, B and / residuals of the common stars of contiguous chips. 
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Chips 


< AV > 


rms 


< AB > 


rms 


< AI > 


rms 


1-2 


0.051 


0.023 


0.043 


0.025 


0.054 


0.047 


1-8 


-0.006 


0.027 


0.033 


0.072 


0.009 


0.020 


2-3 


0.005 


0.036 


0.004 


0.055 


-0.018 


0.031 


2-7 


0.021 


0.070 


0.043 


0.049 


0.014 


0.051 


3-6 


0.046 


0.021 


0.035 


0.036 


0.034 


0.017 


6-7 


-0.032 


0.029 


-0.071 


0.055 


-0.023 


0.026 


7-8 


-0.065 


0.039 


-0.015 


0.043 


-0.043 


0.027 



Table 4. External photometric errors estimated from the artificial star experiments. 



Chip 


V 


Cext{V) 


<Text(B - V) 


a ext (V - I) 


Chip 


V 


Oext{V) 


<Text(B - V) 


<Jext(V — I) 




16.721 


0.003 


0.005 


0.004 


5 


16.873 


0.003 


0.006 


0.003 




17.617 


0.006 


0.009 


0.007 


5 


18.024 


0.007 


0.013 


0.007 




18.663 


0.010 


0.021 


0.013 


5 


19.128 


0.016 


0.037 


0.018 




19.746 


0.027 


0.047 


0.032 


5 


20.111 


0.033 


0.062 


0.043 




20.878 


0.065 


0.150 


0.083 


5 


21.456 


0.113 


0.213 


0.105 


2 


16.738 


0.003 


0.006 


0.004 


6 


16.911 


0.003 


0.006 


0.004 


2 


17.916 


0.006 


0.009 


0.004 


6 


18.049 


0.006 


0.013 


0.006 


2 


18.856 


0.015 


0.021 


0.010 


6 


19.051 


0.015 


0.034 


0.018 


2 


19.989 


0.028 


0.062 


0.027 


6 


20.079 


0.040 


0.085 


0.043 


2 


21.185 


0.094 


0.188 


0.093 


6 


21.459 


0.110 


0.254 


0.116 


3 


16.755 


0.003 


0.006 


0.003 


7 


16.995 


0.003 


0.007 


0.003 


3 


18.064 


0.007 


0.015 


0.009 


7 


18.084 


0.009 


0.018 


0.010 


3 


19.313 


0.019 


0.033 


0.019 


7 


19.144 


0.016 


0.037 


0.013 


3 


20.321 


0.034 


0.086 


0.032 


7 


20.155 


0.035 


0.070 


0.039 


3 


21.318 


0.093 


0.227 


0.126 


7 


21.441 


0.120 


0.246 


0.095 


4 


16.980 


0.003 


0.004 


0.003 


8 


16.804 


0.003 


0.005 


0.003 


4 


17.910 


0.006 


0.013 


0.007 


8 


17.897 


0.007 


0.013 


0.007 


4 


18.942 


0.016 


0.034 


0.012 


8 


18.983 


0.013 


0.036 


0.013 


4 


19.987 


0.031 


0.087 


0.033 


8 


20.088 


0.033 


0.072 


0.038 


4 


21.277 


0.074 


0.231 


0.065 


8 


21.283 


0.095 


0.170 


0.081 



and B ~ 22. The data for the recovered artificial stars of each 
chip have been sorted by the observed magnitude and divided 
into 14 group of about 80 stars. For each group, the median 
observed magnitude V and th e exte rnal errors have been com- 
puted as in lStetson & Harris ( 1988). In Table|4] we report the 
estimated errors for V > 16.5. The values for V < 16.5 are 
less than 0.001 mag. 

2.2. Astrometry 

In order to determine the astrometric solution we have consid- 
ered as reference catalogue the Guide Star Catalogue, Version 
2.2.01 (GSC2.2 STScI, 2001). We have used the Aladin Sky 
Atlas to plot this catalogue on the finding chart of the open 
cluster NGC 3960 field. For each chip we have chosen three 
stars for which we have both the celestial coordinates of the 
GSC2.2 catalogue and the pixel coordinates. We have used 
these stars as reference for the IRAF task ccxymatch; im- 
posing a matching tolerance of 0"3, we have matched a total of 
4224 stars. The IRAF task ccmap has been used to fit a trans- 
formation between pixel coordinates and celestial coordinates 



after applying a 3 a clipping to the data. In our case, a tangent 
plate projection with distortion polynomials has been chosen 
for the sky projection geometry. The scatter plot and the distri- 
butions of the RA and Dec residuals of the transformation are 
shown in Figure |2] for the chips 2 and 3 where the cluster is 
mainly concentrated. We note that the residual distribution is 
highly concentrated with a Gaussian shape in both directions. 
The mean offsets and the rms of the distributions, given on each 
panel, show that the final accuracy is always better than 0"2. 
In order to investigate eventual systematic dependence from the 
position, we have plotted the RA and Dec residuals as a func- 
tion of RA and Dec. The results for the chips 2 and 3 are shown 
in Figure[3]where no noticeable systematic effects are evident. 
Similar results have been found for the other chips. 

3. The Color-Magnitude Diagrams 

The color-magnitude diagrams of all the stars measured in 
our field of view have been obta ined using all sources with 
SHARP parameter JStetsonl[l987l) between -0.8 and 0.8. This 
selection allows us to reject non-stellar objects such as semi- 
resolved galaxies or blended double stars or cosmic rays. The 
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Fig. 2. Scatter plot of the RA and Dec residuals of the transfor- 
mation and distributions of ARA and ADec obtained from the 
comparison of our astrometric solution and the GSC2.2 refer- 
ence catalogue. The vertical and horizontal dashed lines mark 
the mean residuals in RA and Dec. The results are shown only 
for the chips 2 (left panel) and 3 (right panel), where the cluster 
is mainly concentrated. The number of matched stars, the mean 
values and the 1 a rms in arcsec are given on the figures. 

V vs. B — V and the V vs. V — I color-magnitude dia- 
grams of the 39411 selected stars are shown in Figure E|(top 
panels). Due to the Galactic position of the cluster, the dia- 
grams are highly contaminated by foreground and background 
stars making it difficult to discriminate the cluster main se- 
quence. To facilitate the interpretation of the color-magnitude 
diagrams, we have considered the distribution of the stars as a 
function of distance from the cluster center. For this we have 
first calculated the cluster centroid RA ccn = ll' l 50 m 4F;8 and 
Dec con = -55°40'37"4 (J2000), as the median value of the 
celestial positions of stars brighter than V=15.5. Therefore we 
have performed stars counts by using our data to obtain a first 



Fig. 3. Positional residuals, computed as in Figure|2] as a func- 
tion of RA and Dec in the GSC2.2 reference catalogue. The 
results refer to chips 2 (left panel) and 3 (right panel), where 
the cluster is mainly concentrated. 

order approximation of the cluster size. We have derived the 
surface stellar density by performing star counts in concen- 
tric rings around the cluster centroid and then dividing them 
by their respective areas. Figure [5] shows the resulting density 
profile and the corresponding Poisson errors bars. We note a 
flattening of the profile outside 7 arcmin. This density variation 
suggests that most of the cluster members are located within 7 
arcmin, whereas outside the Galactic disk population is dom- 
inant. However, we consider this radius only an approximate 
value where the cluster dominates the field. In fact, as expected 
from dynamical evolution and mass segregation effects, faint 
cluster stars can be located out of this radius. We note, how- 
ever, that our cluster size estimate is significantly greater than 

the value 2. 75 arcmin given for the radius of this cluster by 
I 1 J 1 

Janes ( 198 1 ), where only stars brighter than V = 17 were con- 
sidered. From here on we refer to the above mentioned region 
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inside the circle of 7 arcmin as the "cluster region" and to the 
remaining part of our field of view as the "field region". 

The color-magnitude diagrams obtained using all the stars 
within the cluster region are shown in Figure|4](bottom panels). 
Both the diagrams exhibit the cluster main sequence down to 
V = 1 7 but the faint star main sequence is still hidden by con- 
taminating stars. The Turn Off (TO) point is roughly located at 
(V = 14, B — V = 0.3), while the group of bright stars in the 
red part of the CMD are clearly Red Giant Branch (RGB) stars 
of NGC 3960. 



45 - 
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Fig. 5. Star counts as a function of the distance from the cen- 
troid for all stars detected in the field of NGC 3960. 



3. 1. The Differential Reddening Map 

The cluster main sequence in Figure |4] appears rather broad in 
both diagrams, with higher dispersion in the V vs. V — I CMD 
(lower-right panel). Artificial-star tests indicate that our typ- 
ical photometric errors for V < 17 are crm-v) ~ 0.007 and 
<j(y-i) ~ 0.004; therefore, after ruling out instrumental causes 
and/or photometric errors as possible origin for this effect, we 
have investigated a possible dependence of such broadening 
from the positions of the stars in our field of view. In Fig [6] we 
show the V — I vs V color-magnitude diagrams of 16 subre- 
gions of about 8'6 x 8'.3 in the NGC 3960 field. Celestial coordi- 
nates corresponding to each subregion are given on the top and 
right axes, respectively. The solid line on the color-magnitude 
diagrams is the isochrone of 1.1 Gyr, m etallicity Z = 0.01 
computed by the lPietrinferni et afl (f2003) (see Section^. 

We note that the cluster main sequence is dominant in the 
central regions, while the number of stars in the upper main 
sequence decreases in the external regions. Nevertheless, a dif- 
ferent appearance of the main sequence as a function of the po- 



sition is clearly evident as regards the theoretical model. This 
effect could be due to the fact that the star cluster is present only 
in one small well defined region; such a difference is however 
also evident looking at the blue edge of the field population, 
roughly located from 0.8 to 1.2 in (V — J). We note that this 
shift does not depend from the chip to chip photometric zero 
point offsets because it is present in different CMDs of stars 
falling in the same chip; in addition differences in these CMD 
are not consistent with photometric zero point residuals of ad- 
jacent chips found in Section |2~T1 Thus, we conclude that the 
dependence of the appearance on the spatial positions of the 
stars is real and can be due to differential reddening. 

This conclusion is consistent with several previous studies 
indicating the presence of dust clouds in the NGC 3960 field. 
For instance, a region of 0.116 square degrees, located at the 
RA=ll' l 49 m 10 s , Dec=-55°31'40"(J2000) and (/ =294°15, 
b =6°28), and falling in the north-west part of our field of view, 
has been included in a catal ogue of dark nebulae compiled by 
iFeitzinger & Stuewel |[l984). In this study, many fields of the 
ESO/SRC southern sky survey were examined looking for re- 
gions of the sky where the apparent surface density is reduced 
compared to the surrounding regions. The dark nebula falling in 
our field of view was classified as a cloud with some structure 
and relatively well defined edges. Another dark cloud, falling 
als o in the north - west p art of our field of view, was identified 
bv lHartlev etaD dl986l) using the ESO/SERC Southern J sur- 
vey. The catalogue obtained in this work includes a dark nebula 
of minimum density of size 5' x 3' centered on RA=1 1 M9 m 03 s 
and Dec=-55°40'22"(J2000) (/ =294°16, b =6°13). The dis- 
tances of these dark clouds are not known. 

More recently, a value for the reddening of these re- 
gions was given in the unified catalogue of dust clouds by 
[Dutra & B ica (2002). The reddening values were extracted 
from th e all-sky reddening E(B — V)fir map o flSchlegel et alJ 
(1998), based on a 100 fim dust thermal emission map with 
resolution of ~ 6'. For each region, the authors compare the 
reddening in the nebula direction (E(B — V) ccn ) with the back- 
ground reddening (E(B — V)bck) computed from the average 
of the measured values in four background surrounding posi- 
tions, used as reference. In our field they found E(B — V) ccn = 
0.51 and E(B - V^bck = 0.38, for the first dark nebula and 
E(B - V) celi = 0.59 and E(B - V) hck = 0.55, for the 
second one. These values are significantly different fro m the 
cluster reddening value E(B — V) = 0.29 estimated bv l.Tanesl 
( 198 1 ) as the average of reddening measurements ranging from 
E(B -V) = 0.20 and E(B - V) = 0.34 obtained for 6 giants 
in NGC 3960. Such reddening estimates indicate that strong 
gradients of interstellar absorption are present in our field, con- 
firming our conclusion on the presence of differential redden- 
ing deduced by looking at the color-magnitude diagrams. 

In order to quantify such effects on the color-magnitude 
diagrams we have created two reddening maps, one for the stars 
located at the cluster distance, using the cluster main sequence, 
and a reddening map for the field star population, affected by 
a further reddening spread due to distance. Both maps are used 
to apply an average correction for differential reddening to the 
photometric data. 
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Fig. 4. The top panels show the V vs. B — V and the V vs. V — I color-magnitude diagrams of all the stars measured in our field 
of view, while the bottom panels show the same diagrams of all the stars within the "cluster region". 



As in IPiotto et al.l ifT999). the main sequence reddening map 
has been calculated as follows: we divided our field in 20 x 
20 subregions of about 1(7x1(7 where we have obtained the 
CMDs. Such subregion size is a good compromise to have 
a sufficient number of stars to identify the cluster main se- 
quence, at least within 7 arcmin from the cluster centroid, and 



high spatial resolution to map the differential reddening. We 
have selected as fiducial region the subregion with coordinates 
RA =[ll' l 50 m 28*3, ll ,l 50 m 45 ;i .6] and Dec =[-55°41'13"2, 
— 55°38'50"6] (J2000) where the cluster main sequence is well 
defined and where our estimated cluster centroid is located; we 
have fitted with a spline the points of the CMD where the star 
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Fig. 6. V vs. V — I color-magnitude diagrams of 16 subregions of about 8'6x8'3 in the field of NGC 3 960. T he solid line is a 
representative theoretical isochrone of 1.1 Gyr and metallicity Z — 0.01, calculated by Pietrinferni et al. (2003). 



density is larger, defining the cluster main sequence; for each 
star with 12 < V < 18.2 and 0.4 < (V — I) < 2.4, we have 
calculated the distance from the fitted sequence along the red- 
dening vector defined by the relation Ay = 3.1/1.25 x E(V — 
D, w here E(V - I) = 1.25 x E(B ~ V) JMunari & Carrarol 
Il996h . The relative reddening has been calculated as the me- 
dian value of the V — I component of these distances. 

As in Ivon Braun & Mateol J200ll) . we have obtained the 
reddening map at the cluster distance shown in Figure0 where 
the darker regions correspond to the higher relative reddening; 
for reference, the position of the fiducial region and the cluster 
region are also indicated by the small box and the circle, re- 
spectively. The reddening values in milli-magnitudes, relative 



to the fiducial sequence and corresponding to each subregion, 
are the top number in each pixel of the grid shown in Figure[8] 
The bottom number in each pixel is the number of stars used to 
determine the relative reddening in that pixel. The thicker box 
indicates the position of the fiducial region. We note that higher 
reddening regions correspond to lower star numbers with re- 
spect to the average star number in the surrounding regions. 
Nevertheless, we are aware that the relative reddening values 
are significant for the central regions, where the cluster popula- 
tion is dominant, but are not significant for the external regions 
of our field of view where the stellar population is dominated 
by field stars; in fact, these stars lay at different distances and 
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Fig. 7. Grayscale reddening map of the NGC 3960 cluster members. In this figure, the darker regions correspond to higher relative 
reddening, calculated with respect to the fiducial cluster main sequence. The position of the fiducial region and the cluster region 
are also indicated as the small box and the circle, respectively. 

the fiducial main sequence used to determine the reddening is 
not representative in these regions. 

In order to estimate the differential reddening affecting the 
field stars, we have calculated the separation of each star with 
18.5 < V < 22.0 and 0.0 < (V - I) < 2.8, from the blue edge 
of the CMD in the fiducial region ( von Hip pel et aH l2002). The 
relative reddening affecting these stars has been estimated as 
the mode of the A(V — I) distribution. We note that we are ap- 
proximating the relative reddening to a single value for all the 
field stars, although we are aware that the relative reddening is 
almost certainly spread out in distance. Using a single redden- 
ing value for each subregion is an approximation, since clearly 
we cannot determine the reddening for each star. However, our 
value can be used to determine an average contamination of 
cluster candidates, as will be discussed in Section l5~TI The ob- 
tained results are shown in Figure[9] where the fiducial region 
position is also indicated; the corresponding values in milli- 
magnitudes and the number of the stars used to determine the 
relative reddening are reported in Figure ^3 where the posi- 
tion of the fiducial region is marked by a thicker box. In Figure 
|9j the small circle and the circular line indicate t he po sition 
of the dark nebulae included in the lHartlev et alJ (l986) and 
Feitzinger & Stuewe ( 1984) catalogues. As already mentioned, 
Dutra & Bical J2002h estimated E(B - V) = 0.59 in the re- 
gion corresponding to the small circle. On the other hand, in 
this region we have measured a shift of AE(V — I) = 0.45 



with respect to the CMD blue edge in the fiducial region. Using 
the relation given in lMunari & Carrarol Jl996l) . this shift corre- 
sponds to AE(B — V) — 0.36. By subtracting this value to 
the reddening estimated by Dutra & Bica ( 2002), we find, in 
our fiducial region, E(B — V) = 0.23, that i s cons istent with 
the average value E[B — V) = 0.29 given by Janes ( 198 lj) for 
the central region of the cluster. This external test shows that 
our differential reddening estimate is consistent with previous 
reddening estimates. 

We note that the two reddening maps are qualitatively con- 
sistent with each other, both showing the two regions with high- 
est reddening, located in the north-west and in the south-east 
part of the field, and the lowest reddening region, located in 
the north-east part of the field. In the cluster region the median 
of the differences between the reddening values derived from 
the two methods is Ay = 0.08 mag with a standard devia- 
tion of 0.25 mag. We have used the reddening values relative 
to the cluster main sequence (Figure [8ji to correct magnitude 
and color measurements of the stars within the cluster region, 
where we expect most of the cluster members to be located, 
and the reddening values relative to the blue edge of the CMD 
to correct magnitudes and colors of the stars in the field re- 
gion, where the stellar population is dominated by field stars. 
We find that, within the cluster region, the E(V — I) values 
range from 0.21 up t o 0.78 , assum ing that E(V — I) = 0.36 
(E(B - V) = 0.29 Jjanesl.il 98 llV ) in the fiducial region; cor- 
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Fig. 8. The top number in each pixel of the grid is the reddening values in milli-magnitudes, calculated with respect to the fiducial 
main sequence and corresponding to each subregion, while the bottom number in each subregion is the number of stars used to 
compute the relative reddening in that pixel. The fiducial region is marked by a thick box. 



rections in E(V — I) up to 1 mag have been, instead, applied 
for the most reddened field stars. 

The reddening corrected V vs. V — I CMD in the cluster 
region is shown in Figure^2(^/ f panel). We note that the main 
sequence is better defined with respect to the uncorrected CMD 
shown in the lower-right panel of Figure The Turn Off and 
the RGB stars, which were spread out in Figure^ are now well 
traced. Of course, our procedure computes an average correc- 
tion on each pixel of our map and a residual spread in the CMD 
remains. 

4. Cluster Fundamental Parameters 

Our photometric data together with very recent stellar models 
allow us to determine the cluster parameters, previously esti- 
mated by I.TanesNl98ll) from a sample of only 318 stars with 
V < 16.5. 

In order to find the age and the distance of the cluster, we 
have considered a set of different age theoretical isochrones re- 

^ ^^^^^^^^^ 

cently calculated by Pietrinf erni et alJ ( 120031) . w ith metallicity 
Z = 0.01, as spectroscopically determined by Friel & Janes 
dl993h . By fix i ng the reddening value E(B - V) ~ 0.29, de- 
rived bv ljanesl Jl98ll) . we have vertically shifted the set of red- 
dened isochrones in the V vs. (V — I) CMD of the stars in 
the cluster region. The distance that better fit the upper part of 
main sequence stars (15 < V < 17) corresponds to a distance 
modulus of (V — My)o — 11.35. Since we do not know the in- 



dividual membership we can only constrain the age within the 
range of ages of isochrones limiting most of the stars close to 
the Turn Off. With this criterion the cluster age is between 0.9 
and 1.4 Gyr. 

Our results therefore sugge st an age for NGC 3960 older 
than that given in IJanesI ill 98 ill and a distance of about 1850 
pc, that is slight larger than the value given by the same author 
(d=1660 pc). In the right panel of Figure[^we show the V vs. 
V — I diagram of the stars within the cluster region where the 
adopted theoretical isochrones have been superimposed. 

5. Empirical Cluster Locus and Photometric 
Selection of Candidate Cluster Members 

5.1. Optical Data 

As already mentioned, the reddening corrected V vs. V — I 
CMD computed from the stars within the cluster region is 
highly contaminated by field stars, as shown by the fact that 
only the bright cluster main sequence is clearly visible. In or- 
der to define a sample of candidate cluster members useful to 
study the cluster population, we need to define photometrically 
the complete cluster locus in the CMD. This result is usually 
achieved by superimposing theoretical stellar models to the 
observed CMD where the complete cluster main sequence is 
clearly visible, which is not our case. To check the agreement of 
theoretical isochrones, that suffer from a number of uncertain- 
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Fig. 9. Grayscale reddening map of the field stars in the NGC 3960 field. In this figure, the darker regions correspond to the 
higher relative reddening, calculated with respect to the blue edge of the CMD in the fiducial region. The position of the fiducial 
region is indicated as the small box. The small circle and the circular line indicate the position of the dark nebulae included in 
the Hartley et al. (1986) and Feitzinger & Stuewe (1984) catalogues. 



ties ( von Hip pel et all l2002t iGrocholski & Saraiedinil l2003t> . 
with our data, we have empirically recovered the cluster main 
sequence using a statistical subtraction of the CMD density dis- 
tribu tion, as described belo w. A similar approach has been used 
bv lvon Hippel et alJJ2002l) . 

We have used the field region to perform a statistical sub- 
traction of the contaminating field stars in the CMD derived 
in the cluster region. The result of this subtraction defines 
the locus of the cluster main sequence JChen et all Il998l 
Ivon Hippel etaill2002tlBaume et alll2003h . 

Using our reddening corrected photometric data, we have 
defined a "clean field region " by selecting from the "field re- 
gion" all those subregions in FigurefTolwhose relative redden- 
ing correction was between —0.150 and 0.200, corresponding 
to the range of reddening, estimated for the field stars within 
the "cluster region". We have chosen these limits in order to 
consider only field stars affected by an absorption similar to 
that of the cluster stars. For both the cluster and the field re- 
gions, we have constructed a grid of boxes of AV = 0.125 and 
A(V — I) = 0.05 on the V vs. V — I plane. Thus, we have 
built a greyscale CMD density map of the cluster by subtract- 
ing the unreddened CMD density map of the field stars, from 
the unreddened CMD density map computed for the "cluster 
region". Figure ITU shows the resulting map where only the 
boxes showing a 1 a excess are shown. This excess should fol- 



low the cluster main sequence. The cluster main sequence is 
clearly present down to V ~ 18, while the subtraction is nois- 
ier for V > 18.5, the region of the cluster main sequence more 
strongly contaminated by field stars, where the main sequence 
exhibits a lower star density; as we have already mentioned, our 
relative reddening correction suffers from the uncertainty due 
to the spread of the reddening of individual stars, which pre- 
vents us from obtaining an accurate field star CMD estimate. 
This explains the lower star density in the low mass cluster 
main sequence and the noise of the bluer field stars. However, 
in this subtraction the faint cluster population (V > 18.5) is 
well separated by the field star population this allowing us to 
have a well defined empirical locus of the cluster main se- 
quence. 



In Figure 1 1 2.1 the 0.9 and 1.4 Gyr isochrones 
JPietrinferni et all 120031) have been superimposed onto 
the CMD density map using the cluster parameters found 
in Section |4] The good agreement between the empirically 
recovered cluster main sequence locus and the adopted theo- 
retical model down to at least V ~ 18, confirm both that the 
parameters determined in Section |4] are suitable to the cluster 
and that the adopted theoretical models reproduce the cluster 
main sequence over the mass range for which the IMF will be 
fitted. 
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Fig. 10. The top number in each pixel of the grid is the reddening values in milli-magnitudes, calculated with respect to the blue 
edge of the CMD in the fiducial region and corresponding to each subregion, while the bottom number in each subregion is the 
number of stars used to compute the relative reddening in that pixel. The fiducial region is marked by a thick box. 



We have used both the empirical cluster locus and the men- 
tioned theoretical isochrones to define a photometric cluster 
member sample. To this aim, we have defined a strip in the 
observed CMD; the lower and upper limits take into account 
the spread of the reddening determined above. Furthermore, 
the upper limit of our strip is displaced upward by 0.75 mag 
in order to include binary stars. We have defined as candidate 
cluster members a total of 2119 stars lying in the cluster re- 
gion that, within their photometric errors, belong to this strip. 
The photometric/astrometric catalogue of these candidate clus- 
ter members is given in Table 5 4 , where we report RA and Dec 
(J2000) coordinates in decimal degrees, an identification num- 
ber for each star, the V, B, I magnitudes and their uncertain- 
ties. 

5.2. IR Data 

As already stressed in the previous Sections, our field of view 
is strongly affected by differential reddening, thus making un- 
certain the field star CMD estimate and the study of the clus- 
ter population. To have an external check of our results, we 
have used near-infr ared JHKj, data from 2MASS survey of the 
same field of view JCarpente3.l200ll) . since IR bands are not as 

4 available in the electronic form at the CDS via anony- 
mous ftp to cdsarc . u-strasbg . fr (130.79.128.5) or via 

|http : / /cdsweb .u-strasbg.fr/cgi-binl 



much sensitive to reddening effects as optical bands. To have 
a sample as reliable as possible, only stars detected in all the 
JHKs bands, with magnitudes measured from point spread- 
function fitting or aperture photometry have been selected from 
the 2MASS catalogue. 

The first step has been to construct IR CMDs of the stars 
in the cluster region. Figure^]shows the J vs. J — K$ CMD 
of 1162 stars detected in this area; we can see that the cluster 
main sequence is revealed from higher star density sequence 
down to J ~ 16.5, although it remains rather broad probably 
due to a large scatter in the colors of the low Signal-to Noise 
ratio sources. 

In this diagram, the clump of stars at J ~ 11 and (J — 
K s ) 0.7 are certainly Red Giant Branch stars of NGC 3960, 
while the lower density stars in the almost vertical sequence 
are mainly giant field stars. As already done with optical data, 
we have built the greyscale density map of the cluster with IR 
data by subtracting the field CMD density map from the J vs. 
J — Ks CMD density map computed for the "cluster region". 
In order to estimate an optical-independent field star distribu- 
tion, we have defined as "field region " the farther annulus area 
with radius R between 30' and 45' from the cluster centroid. 
Figure [H] shows the IR resulting map where only the boxes 
having a 1 a excess are shown. We note as the cluster main 
sequence is well defined at least down to J = 15.8, that is the 
2MASS completeness and reliability magnitude limit. In or- 
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Fig. 11. Left panel: The V vs. V—I diagram of the stars within the cluster region, corrected for differential reddening. Right panel: 
The V vs. V — I diagram of the stars within the cluster region, corrected for differential reddening, as compared to the Pietrinferni 
et al. (2003) isochrones of age 0.9 and 1.4 Gyr, for the metallicity Z = 0.01. Using the color excess of E(B-V) = 0.29 iTanesl 
Il98li a distance modulus of (V - M v ) = 11.35 mag has been derived. 



der to verify the cluster parameters, we have superimposed in 
Figur e H4l the 0.9 and 1.4 Gyr isochrones of P ietrinferni et alJ 
(2003), adopted also in the optical case, using the parameters 
that we have determined in Sec tion |4] and the reddening law 
derived by Cardelli etal.ldl989l) . We again find a good agree- 
ment between theoretical models and observational data down 
to J ~ 15.5, while for fainter stars, the model do not repro- 
duce the empirical cluster main sequence; this can be due both 
to the s trong uncertainties affecting IR color-temperature re- 
lations dGrocholski & Saraiedinii 120031) and to the photomet- 
ric errors in this magnitude range. However, the agreement of 
the adopted theoretical model with the bright cluster main se- 
quence confirms again the reliability of the cluster parameters 
found in Section^ 



Once the empirical cluster main sequence has been recov- 
ered from the IR data, we have applied, as done with the opti- 
cal data, a photometric selection of candidate cluster members. 
Using the adopted theoretical isochrone for the upper main se- 
quence and a fiducial line for the lower one, we have defined 
a strip in the CMD, including the boxes about the main se- 
quence with a 1 a excess; in this way we have selected a sam- 
ple of 749 IR cluster candidates within the "cluster region " and 
belonging, according to their photometric errors, to this strip. 
The photometric/astrometric catalogue of these candidate clus- 



ter members is given in Table 6 5 , where we report RA and 
Dec (J2000) coordinates in decimal degrees, an identification 
number for each star, the V, B, I magnitudes and their uncer- 
tainties. 

6. Luminosity and Mass Functions 

One of the main objectives of this study is to derive the stellar 
mass function down to the limiting magnitude reached in this 
survey, that is V ~ 22. 

In Section 15.11 we have found a sample of optical candi- 
date cluster members, but such photometrically selected sam- 
ple contains the contribution of foreground and background 
sources that has to be subtracted in order to estimate the number 
of stars belonging to NGC 3960. Since neither spectroscopic 
observations, nor proper motions are available for this cluster, 
we have used the following statistical approach to estimate the 
number of the cluster members. 

First, we have computed the total distribution of the optical 
candidate members lying in the "cluster region", as a function 
of the V magnitude; second, we have computed the V distribu- 
tion of the candidate members falling in the "clean field region" 
defined in Section lBTl unreddened relatively to the fiducial re- 
gion; then, this field distribution has been normalized to the 

available in electronic form at the CDS via anonymous 
ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 

|http : / /cdsweb . u-strasbg . f r| 
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Fig. 12. Grayscale CMD density map of the cluster obtained by 
subtracting the V vs. V — I CMD density maps, computed for 
the "clean field region" and for the "cluster region". Only boxes 
with alcr excess are shown. The solid and dashed lines are, 
respectively, the 0.9 and 1.4 Gyr isochrones with metallicity 
Z = 0.01, computed by Pietrinferni et al. (2003), scaled to the 
observed CMD using the cluster parameters given in Section|4] 



"cluster region" area and then subtracted to the total one. The 
resulting V luminosity function of the cluster is shown in the 
top panel of Figure ^] where the values of the absolute mag- 
nitude (My), corresponding to the visual magnitudes are also 
indicated in the top axis. We note that for V > 18.5, the distri- 
bution is in some case consistent with the zero value and in one 
case, is even negative(!). This result is consistent with the low 
star density found in Figure ^] in the same magnitude range. 
As already mentioned, this finding is very likely due to the crit- 
ical field star distribution estimate about V = 19.5, where the 
contribution of the field stars is dominant. For V > 20, the 
deficiency of stars is probably due to incompleteness of our 
catalogue that, based on the artificial star tests, is 100% com- 
plete down to V ~ 20. In addition, due to dynamical evolution, 
a percentage of such faint cluster stars is expected to be found 
also outside our estimated cluster region, that is where we have 
estimated the field star distribution, and this makes it difficult 
to have a correct estimate of the cluster luminosity function in 
this magnitude range. 

An independent IR luminosity distribution of the cluster 
has been obtained using the same statistical approach and the 
sample of IR candidate cluster members obtained in Section 
15.21 As for the optical case, the contribution of field stars has 
been computed using the "field region" defined in Section l5~2l 
The resulting J luminosity function is shown in the bottom 
panel of Figure [21 Also in this case, the result is reliable for 
J < 15.8, that is in the magnitude range where the 2MASS 



Fig. 13. J vs. J — K s CMD obtained from the 2MASS cata- 
logue of the stars detected within 7 arcmin from the cluster cen- 
troid. The diagram shows only stars revealed in all the JHKs 
bands with magnitudes measured from point spread-function 
fitting or aperture photometry. 

catalogue is complete. This limiting magnitude corresponds to 
V ~ 17, for a main sequence star at the cluster distance. 

In order to compare these results, obtained from two inde- 
pendent surveys, we have transformed the V and J luminos- 
ity functions into the corresponding mass functions. We have 
used in b oth cases, the Mass-Lum inosity Relation (MLR) com- 
puted by Pietr inferni et alJ J2003). The derived mass functions 
are shown in Figure^] As already stressed, the results are re- 
liable for M > 1M©, where both surveys are complete. For 
this mass range, we have computed a linear fit to both distri- 
butions and have found that the two functions have a slope of 
2.95 ± 0.53 and 2.81 ± 0.84, for the V and J luminosity func- 
tions, respectively. 

The agreement of the mass function obtained from the op- 
tical data with the one obtained from IR data indicates that the 
applied corrections for differential reddening and, therefore, 
the field star distribution estimate, are reliable. Nevertheless, 
the very large errors corresponding to the optical mass function 
in the range M < IMq are the result of the already mentioned 
problems in this mass range that prevent from having a reliable 
estimate of the low mass function of NGC 3960. 

7. Summary and Conclusions 

This Chapter presents photometry and astrometry of the stars 
in the field of the open cluster NGC 3960, falling in the 34 x 33 
arcmin square field of the WFI camera of the MPG/ESO 2.2 m 
Telescope. Our survey reaches a limiting magnitude V ~ 22; 
based on artificial star tests, it is 100% complete down to 
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Fig. 14. Greyscale density map of the cluster obtained by sub- 
tracting the J vs. J — Ks CMD density map, computed for the 
"field region" (cf Section l5~2l . from the J vs. J — K$ CMD 
density map, computed for the "cluster region". Only boxes 
with a 1 (j excess are shown. The solid and dashed lines are, 
respectively, the 0.9 an d 1.4 Gyr isochrones w ith metallicity 
Z = 0.01, computed bv lPietrinferni et alJ J2003I) . scaled to the 
observed CMD using the cluster parameters given in Section|4] 

V ~ 20, where the photometric accuracy is better than 6%. 
Comparison of the positions of the stars in our survey with 
those of the GSC 2.2, used as reference catalogue, indicates that 
the final astrometric accuracy is <0"2. 

The photometric data have been used to construct the V vs. 
B - V and the V vs. V — I CMDs, where the high contribu- 
tion of the disk population allows us to identify only the cluster 
bright star main sequence, that appears rather broadened due to 
strong differential reddening. A detailed spatial study of the V 
vs.V — I CMD allows us to evaluate the relative reddening of 
20 subregions, each of about l'.7xl'.7, with respect to a fidu- 
cial region corresponding to the cluster centroid. The relative 
reddening has been computed from the cluster main sequence, 
for the stars falling in the cluster region, and using the blue 
edge of the CMD for the field stars, which, being at distances 
different from that of the cluster, are affected by different red- 
dening. Our results indicate that, within the cluster region, the 
E(V — I) values range from 0.21 up to 0.78; corrections in 
E(V — I) up to 1 mag have been, instead, applied for the most 
reddened field stars. 

The final reddening map traces regions of low star density, 
that are consistent with two areas, indicated in literature as dark 
nebulae, that are regions of the sky where the apparent surface 
density of stars is reduced compared to surrounding regions. 
Comparison of our reddening estimate in the smallest of these 
dark nebulae, with the one reported in literature, allows us to 



Fig. 15. Luminosity functions of NGC 3960, corrected for the 
field star contribution, obtained from optical data (top panel) 
and from IR data {bottom panel). The dashed lines indicate the 
zero value of the distributions. 



estimate a color excess E(B — V) of about 0.23 in the region 
corresponding to the clust er centroid. T his value is very similar 
to the value 0.29 given bv LfanesI dl98lh for the cluster. 

Using the reddening corrected CMD and recent stellar 
models, we have found that the cluster has an age between 0.9 
and 1.4 Gyr, its distance modulus is (V — Mv)o = 11-35 and 
it is located about 1850 pc from the Sun. 

In order to minimize the field star contamination in the 
CMD, a statistical subtraction of the CMD density map has 
been performed to empirically recover the whole cluster main 
sequence. Using the resulting density map, we have been 
able to compare the empirical cluster main sequence with the 
adopted theoretical models, thus allowing us to verify the re- 
liability of our estimated cluster parameters. Furthermore, this 
empirical cluster main sequence has been used to identify a 
strip in the optical CMD where cluster stars are expected to be 
located; 2119 candidate cluster members fall in this strip. 

In order to have an external check of our results, we have 
built an analogous CMD density map using IR data from the 
2MASS catalogue. Taking advantage of the availability of the 
entire sky coverage of the 2MASS catalogue and in order to 
avoid the above mentioned dark nebulae, a different control re- 
gion for field stars, located farther than the one used in the op- 
tical case, has been used. The obtained empirical cluster main 
sequence in the J vs. J — Ks plane has been compared to the- 
oretical models using the cluster parameters estimated in this 
work. The result of this comparison shows, again, the correct- 
ness of the estimated cluster parameters. As in the optical case, 
the IR empirical main sequence has been used to define a strip 
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Fig. 16. Comparison of the NGC 3960 mass function obtained 
from the optical data (filled points) and from IR data (empty 
triangles). The solid line is the power law fit to the mass func- 
tion obtained from the V luminosity function, while the dashed 
line is the power law fit to the mass function obtained from the 
J luminosity function. £(M) values are given in number per 
logarithmic mass unit. The slopes of the power law fitting with 
the rms residuals of the least squares fit, as obtained from the 
two distributions, are also indicated. 



in the J vs. J — Ks CMD and therefore a sample of IR candi- 
date cluster members. 

From the two samples of candidate cluster members, ob- 
tained from the optical and IR surveys, and using the respec- 
tive control field, the luminosity distributions for the cluster, 
as a function of the V and J magnitudes, have been com- 
puted. Using the same MLR, the two distributions have been 
transformed into two independent mass function determina- 
tions of the cluster. In the mass range where the two data set 
are complete (M < 1M Q ), the two mass functions have been 
fitted to a power law having indices ay = 2.95 ± 0.53 and 
a j = 2. 81 ±0.84 in V and in J, respectively, while the Salpeter 
mass function in this notation has index a = 2.35. Our values 
are bo th consistent wi th the data presented in the a vs. log M 
plot of Kroupal J2002h for other open clusters. 

The good agreement of the two distributions, furthermore, 
suggests both that the reddening correction of our catalogue 
and the estimates for the field star contamination are reliable 
and correct. However, the low mass range is affected both by 
completeness problems as well as by likely mass segregation 
effects preventing a reliable derivation of the cluster luminosity 
and mass functions for M < 1M Q . 

A future spectroscopic study of this region is planned in 
order to measure the radial velocities crucial to have an inde- 



pendent membership criterion and to determine individual red- 
dening of each star. 
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